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ABSTRACT 

We introduce a simple linear equation relating the line-of-sight peculiar velocity and den- 
sity contrast correlation functions. The relation, which we call the Gaussian cell two-point 
"energy-like " equation, is valid in the distant-observer-limit and requires Gaussian smoothed 
fields. In the variance case, i.e., at zero lag, the equation is similar in its mathematical form 
to the Layzer-Irvine cosmic energy equation. (3 estimation with this equation from the PSCz 
redshift galaxy survey and the SEcat catalogue of peculiar velocities is carried out, returning 
a value of (3 = 0.44 ± 0.08. The applicability of the method for the 6dF galaxy redshift and 
peculiar motions survey is demonstrated with mock data where it is shown that beta could 
be determined with « 5% accuracy. The prospects for constraining the dark energy equation 
of state with this method from the kinematic and thermal Sunyaev-Zel'dovich cluster surveys 
are discussed. The equation is also used to construct a nonparametric mass density power 
spectrum estimator from peculiar velocity data. 

Key words: galaxies: clusters: general - galaxies:distances and redshifts - cosmology: theory 
- large-scale structure of Universe - cosmological parameters 



1 INTRODUCTION 

In the linear regime of the gravitational instability scenario the un- 
derlying mass distribution is directly traced by the galaxy peculiar 
velocities. Measurement of the radial component of the galaxy pe- 
culiar velocities, the only component that one can easily observe, is 
carried out with one of many available techniques, the most com- 
mon among which are the Tully-Fisher-like methods (e.g., Tully 
& Fisher 1977, Faber & lackson 1976). Normally, these meth- 
ods exploit a well defined intrinsic relation between two or more 
of the galaxy 1 observed properties that facilitates establishing its 
actual distance from the observer. The estimated distance is then 
used together with the measured galaxy redshift, to determine the 
galaxy radial peculiar velocity. Assuming an irrotational flow on 
large scales and the knowledge of Q m (the cosmological mass den- 
sity parameter), it is straightforward to use the measured radial pe- 
culiar velocities to recover the full underlying mass overdensity 
(Bertschinger & Dekel 1989, Dekel, Bertschinger & Faber 1990, 
Zaroubi 2002, Zaroubi et al. 2002). In addition, the same mass- 
density could be probed by galaxy redshift catalogues assuming a 
simple linear biasing scheme that connects it to the spatial galaxy 
distribution (Kaiser 1984, Bardeen et al. 1986). To date, galaxy red- 
shift and peculiar motion surveys are the main tools with which as- 
tronomers explore the distribution of matter in the nearby universe. 



1 It should be noted that some methods are not based on galaxy properties 
and use other "extra-galactic objects", e.g., supernovae-Ia. 



Since the two types of data, galaxy peculiar velocity cata- 
logues and galaxy redshift surveys, probe the underlying mass dis- 
tribution comparing the two provides a simple and powerful test 
on the paradigm of gravitational instability and gives a model in- 
dependent measurement of [3, the ratio between the linear growth 
factor, fiyim) (~ S^m 6 ) an d the linear biasing factor of the galaxy 
population. In most cases the comparison is either performed by de- 
riving galaxy peculiar velocities from the galaxy density field and 
confront them with the measured velocities point-by-point, an ap- 
proach usually called "velocity-velocity" comparison (e.g., Davis, 
Nusser & Willick 1996, Willick & Strauss 1998, Zaroubi 2002, 
Zaroubi et al. 2002). Or by adopting the so called "density-density" 
approach in which the velocity data is used to infer the full mass 
density field (Bertschinger & Dekel 1989, Zaroubi 2002) and com- 
pare it with the galaxy distribution (e.g., Sigad et al. 1998, Zaroubi 
et al. 2002). With the exception of the POTENT algorithm (see e.g., 
Sigad et al. 1998), all the comparison methods yield a low value of 
(3 consistent with fi m « 0.3 and bias factor of w 1. 

Another approach to the comparison that doesn't fit into 
the two general classes outlined earlier, is the one proposed by 
Juszkiewicz et al. ( 1 999) who start from the pair conservation equa- 
tion (Peebles 1980) and evolve it further to the quasi-linear regime 
of gravitational instability. In the pair conservation approach, which 
yields a value of f3 that is consistent with the one derived from 
velocity- velocity analyses (Ferreira et al. 1999, Feldman et al. 
2003), a relation between the mean pairwise velocity at a certain 
separation and the density correlation function is derived. The com- 
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parison in this approach is significantly simplified by avoiding the 
spatial point-by-point matching required in previous methods thus 
reducing the noise involved. The Juszkiewicz etal. (1999) approach 
is similar to the one developed in this paper except that we are inter- 
ested in the variance of the peculiar velocity at a given smoothing 
scale rather than the pairwise peculiar velocity at a given separa- 
tion. 

In this study, we derive a very simple, model independent and 
linear relation, valid in the distant observer limit, between the over- 
density and peculiar velocity two point correlation functions as- 
suming Gaussian smoothing. The method is first used to construct a 
nonparametric estimator of the mass-density power spectrum. Then 
the paper concentrates on the relation between the variance (2-point 
correlation at zero distance) of the two fields. This relation basically 
reduces the comparison between the catalogues to two numbers al- 
lowing a robust extraction of the parameter 0. The proposed equa- 
tion is especially suited to future peculiar velocity data sets like, 
the 6dF which will measure the the peculiar velocities of 15, 000 
galaxies with their D n ~ a relation up to 150 Mpc/h distance. 

Currently, the main sources of error in the redshift-peculiar 
motion comparison are the large random and systematic uncer- 
tainties carried by the peculiar motion measurements, for exam- 
ple the Tully-Fisher-like relations has an inherent uncertainty of 
~ 15 — 20% of the distance. Data obtained with more accurate 
distance indicators do exist (Tonry 1991, Riess, Press & Kirsh- 
ner 1995), unfortunately however, either they are not at significant 
distances, e.g., the Surface-Brightness-fluctuations method, or they 
reach large distances but are too sparse, e.g., the Supernovae-Ia 
data. In the future, by using the thermal and kinematic Sunyaev- 
Zeldovich (SZ) effect (Sunyaev & Zeldovich 1972), both the accu- 
racy of the peculiar velocity measurement and the spatial coverage 
of the data are expected to increase dramatically where the uncer- 
tainty is expected to amount to an absolute error of « 150 km s -1 
(e.g.,Diaferio et al. 1994) and the number of observed objects to 
reach ~ 10 4 clusters. The main difficulty in this case, assuming a 
reasonable control over the systematics, will be posed by the large 
mean separation between the galaxy-clusters observed with the SZ 
effect. 

The power of the method proposed here is that it reduces the 
contribution of the measurement noise and sparseness of the sam- 
ple to a bare minimum. The paper is organized as follows: section|2| 
presents the main theoretical formulae. In section|5|the method is 
applied to the PSCz redshift galaxy catalogue and the SEcat pecu- 
liar velocity survey. In section|4]the applicability of the method to 
future surveys, e.g., the 6dF galaxy survey and the kinematic and 
thermal SZ cluster survey, is discussed. 



2.1 The basic relation 

Consider a radial peculiar velocity field vi oa (r) measured in a very 
distant patch of the sky smoothed with a Gaussian window func- 



tion with scale R a , Wr s (t) = (2ttR 3 



-3/2 



exp(-r727?f). In 



the limit of R c -C R, where R c is the correlation radius of peculiar 
velocities and R is the distance of the patch from the observer. A 
smoothed radial field within a given observed volume can be writ- 
ten as, 
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Wfl s (fc)e X p(-ik-x)d 3 fc, (1) 



where the superscript S refers to values smoothed with a Gaussian 
kernel of radius R s , f ; os is a unit vector along the line-of-sight, and 
Wr s (k) is Fourier transform of the smoothing kernel. 

The two-point correlation function of the Gaussian smoothed 
line-of-sight galaxy peculiar velocity is: 



(Cr(r;R s )) 
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Where f\. is the mass density power spectrum and r is the radius 
vector separating between any two points. Since there are two in- 
dependent directions that appear in eq.|5|one can't invoke symme- 
try arguments in order to proceed. However, since in the distant- 
observer-limit the line-of-sight direction is approximately constant 
across the observed volume and independent of the direction of r, 
one can average over all possible directions of r relative to k (arbi- 
trary direction) by integrating equation [3] with | J dfi, where fi 
is cosine the angle between the two vectors r and k: 
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Where () M is an average over statistical ensemble and over \i, 
jo(kr) is the zero order Spherical Bessel function and r = |r| 

Assuming statistical isotropy for the velocity field, i.e., sym- 
metry between the line of sight and the other two orthogonal direc- 
tions one obtains the following equation: 



<e s (r;i? s )> M = 



p 2 H 2 



3(2^ 



. , -£Wg(R s )jo(kr)d d k. 



(5) 



With the factor 3 coming from the symmetry argument. Now to 
the last step in the calculation, for a Gaussian smoothing kernel, 
i.e., Wr s (k) = exp( — k 2 R 2 /2), the derivative of the line of sight 
velocity two point correlation function with respect to R s , yields: 
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2 THEORETICAL DERIVATIONS 

In this section, we first derive the main theoretical relation (sub- 
section [23 and show how it could be used to estimate the matter 
power spectrum from peculiar velocity data (subsection l2.2t . Then 
in subsection l2.3.1l the variance component of the main relation is 
used in order to estimate the value of f3 from comparison between 
galaxy redshift surveys data and peculiar velocity data. The j3 mea- 
surement error for a typical case is derived in subsection 12.3.21 
The derivation is performed within the framework of linear gravi- 
tational instability, under the assumption of statistical homogeneity 
and isotropy. 



Here £(r;R s ) is the two point correlation function of the 
smoothed densities. Notice that eq. Qcan only be obtained when 
dWR s (k) 1 'dR a oc k 2 Wn s (k), strictly valid only with Gaussian 
smoothing kernel 2 . 

Obviously, for a given 3 dimensional peculiar velocity field the 
two point correlation function of the Gaussian smoothed velocity, 



v , is related to its density counterpart through the equation, 



dR a 



= (v s (x)-v s (x- 



(8) 



2 There are other functions that satisfy this relation but they do not satisfy 
the requirements of smoothing kernels 
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which is similar to equation [7] without the factor of 3 and with no 
need for averaging over /i. 

It might be easier to interpret equation Q in its integral form, 
where the integral is performed over the smoothing radius. Let Ri 
and R2 the two smoothing radii that bound our integral, therefore, 

H2 p( r - ra d 3 R„ 

(10) 
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In the r = limit, the left-hand side of equation fTtjl describes the 
mean change in the kinetic energy associated with the smoothed 
velocity due to the variation of the smoothing radius, whereas the 
right-hand side depicts the 3-dimensional integral of the density 
variance of the smoothed field over the smoothing scale. The right- 
hand side term is very similar to the normal potential energy ex- 
cept that R s does not represent a proper distance between points. 
In other words, the variation in the kinetic-like energy comes from 
the "potential-energy-like" behavior of the modes corresponding to 
the scales between Ri and R2, with exponentially decreasing con- 
tributions from larger and smaller scale modes. If r 7^ then the 
interpretation is not as simple but still the left- and right-hand sides 
correspond to a sort of a two-point kinetic-energy and potential- 
energy, respectively, with the same scales contributing to the mod- 
ification as before. Therefore, equation fTol is an "energy equation" 
of sorts as it describes the two point "potential-like" and "kinetic- 
like" energy partition within a Gaussian window function. Dimen- 
sional arguments significantly restrict the mathematical form of the 
relation between the velocity and the density 2 points correlation 
functions. Therefore, it is not surprising that the theoretical relation 
shown by equation Q is very similar to the Irvine-Layzer cosmic 
energy equation which describes how the energy of the Universe 
is partitioned between kinetic and potential energy (Irvine 1961, 
1965; Layzer 1963, 1966; see also Mo, Jing & Borner 1997). 

2.2 The mass density power spectrum 

The main approach currently used to measure the mass-density 
power spectrum from peculiar velocity data is the likelihood 
method introduced by Zaroubi et al. (1997; see also, Freudling 
et al. 1999 and Zaroubi et al. 2001) in which a theoretical power 
spectrum with few free parameters and a noise model are assumed. 
Since in this method the data is forced to fit a specific power spec- 
trum shape, an inaccurate description of the noise model could 
propagate to large scales and contaminate the measured power 
spectrum. Indeed the results obtained with this method have been 
yielding unrealistically high amplitude of the mass-density fluctu- 
ations power spectrum (consistent with Q m > 0.6). Therefore, di- 
rect nonparametric methods for power-spectrum estimation from 
peculiar velocity data are needed. 

The question we pose here is: can equation [7] be used to di- 
rectly estimate the mass power spectrum from peculiar velocity 
data? In the ideal case in which the uncertainties in the measure- 
ment are neglected and the data extends to infinity and samples the 
universe very accurately, the answer is clearly yes. 

A good point from which to start the derivation of the power 
spectrum estimator is the relation between the power spectrum and 
the smoothed density two point correlation function, 
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r from zero to 00 after multiplying with jo(kr)r 2 , one obtains the 
following relation, 

d(er(r;R s )) 



67T 

H^Rs 



dR 3 



-jo{kr)r dr 



fP k e~ k R °. (12) 



Where the orthonormality of the spherical Bessel function is used 
(e.g., Arfken & Weber 2002). The left hand side of eq. H2l is a quan- 
tity that can be directly measured from the velocity data; whereas, 
the right-hand side shows the estimated quantity. Since we use 
smoothed velocity data, the power spectrum is determined up to 
a factor of f 2 (Q, m ) and with a resolution that cannot exceed the 
scale imposed by the Gaussian smoothing. 

For a real application, the discrete and noisy nature of the 
data should be taken into account, i.e., some of the steps leading 
to eq. H2l has to be modified. For example, to maintain the orthogo- 
nality of the spherical Bessel functions on a finite spherical volume 
one has to impose appropriate boundary conditions (see Fisher et 
al. 1995 & Zaroubi et al. 1995 for examples). However, the main 
hurdle for this direct approach to power spectrum estimation is the 
noise contribution; this issue is deferred to a future work. 

2.3 Estimation of f3 

2.3.1 Estimator 

Eq.0is the most general relation derived in this paper. However, 
in order to use it to estimate the value of f3, it is simpler to restrict 
ourselves to the relation between the density and velocity variances, 
namely, apply the equation in the limiting case of r = to yield: 



-~/3 2 H$R s o-j(R s 



(13) 



do 2 (R s ) 
dR 3 

Where a 2 and 0$ are the peculiar velocity and density-contrast 
variances, respectively. 

The numerical calculation of da 2 (R s ) /dR s is straightfor- 
ward. One has to smooth the measured velocity field with a Gaus- 
sian window, calculate its variance and obtain its derivative by finite 
differencing (see subsection l2. 3 . 2l for a similar explicit calculation). 
The right-hand-side of equation fT3l is obtained from the galaxy red- 
shift catalogue by taking the variance of the smoothed real-space 
density field. 

The proposed estimator requires no heavy data manipulation 
and is easy to calculate. Due to the smoothing involved, the estima- 
tor is robust with regard to instabilities caused by the large random 
noise. In addition, to avoid the cosmic variance contribution to the 
error analysis, the comparison between the two types of data sets is 
performed within the same region of space. Both features, simplic- 
ity and stability, render the estimator very appealing to use. 



2.3.2 Noise 

The contribution of the measurement error to the estimator in equa- 
tion^] is readily calculated with the following discrete approach. 
Let e(r-i) be the noise associated with particle i, then the smoothed 
noise is, 



(r0=> XTi)Wn m (r i -T l ) 



1 



(14) 



Substituting this into equation and integrating over the variable 



Subsequently, the expectation value of the noise two point cor- 
relation is, 

{e s {v,)e s {v 3 )) = J2{t 2 (n))W Rs {n - vAWr^tj - r,). (15) 

1 
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The last equation assumes that the measured errors are statistically 
uncorrelated. 

We now require that n = ij and sum over all the data points. 
The expectation value of the noise contribution to the variance of 
the velocity is: 

a 2 N (R s ) = l^( e 2 ( ri )}VKi( ri - ri ). (16) 

i,l 

Therefore, the noise variance that adds to the left hand side of equa- 
tion ll3l is readily obtained by finite differencing: 

da 2 N {R s ) _ a 2 N (R s + AR S ) - a 2 N {R s ) 

dR s ~ AR a ■ ( ' 

The contribution of the noise variance to the right hand side of 
eq. is typically small and is neglected here. However, it is 
straightforward to account for in the case of unusually noisy data. 



3 COMPARISON BETWEEN THE SECAT AND THE 
PSC2 CATALOGUES 

In this section equation^|is employed for comparison between the 
PSCz galaxy redshift catalogue (Saunders et al. 2000, Branchini et 
al. 2000) and the SEcat galaxy peculiar velocity catalogue (Zaroubi 
et al. 2002) which is a combination of the two homogeneous pecu- 
liar velocity catalogues, the SFI catalogue of spiral galaxies (Gio- 
vanelli et al. 1998, Haynes et al. 1999) and the ENEAR catalogue 
of early-type galaxies (da Costa et al. 2000). The SEcat catalogue 
extends to a distance of about 70 /i Mpc and the PSCz goes to 
about twice of that. Therefore, in order to avoid cosmic variance 
contamination of the measurement the comparison between the two 
is restricted to the closer distance. 

Prior to applying the method to the actual data, however, one 
needs to address the question of whether it is realistic to expect a 
reliable estimation of the value of with noisy and close by cata- 
logue such as the SEcat. Hence, the next subsection is dedicated to 
testing with mock catalogues how robust our estimator is. 

3.1 Testing with mock data 

The density and peculiar velocity mock catalogues used in this sec- 
tion are derived from the 3.2 7i Mpc resolution reconstruction of 
the density field from the PSCz galaxy redshift catalogue (Bran- 
chini et al. 2000), where the peculiar velocity field is obtained us- 
ing linear theory from the galaxy redshift space positions assuming 
a value of /3 — 0.5. The mock SEcat peculiar velocity catalogue 
has the same distances and number of points the real SEcat has, but 
with the velocities of the PSCz reconstructed velocity field. Ob- 
viously, it would have been better to use a full nonlinear N-body 
simulation with which to test the method. However, since the posi- 
tions of the actual measured velocities are controlled by the specific 
distribution of the galaxies in the nearby universe we choose to test 
the method with data that has the same spatial distribution as the 
real universe, albeit the lack of full nonlinearity. Given the heavy 
smoothing involved in the analysis this way of assigning veloci- 
ties to mock data is satisfactory - for testing the method with full 
nonlinear simulation see section|4| After assigning the velocities to 
the noise free mock data, we generate 30 mock catalogue with the 
random errors added to their distance and velocity values in con- 
cordance with the observational uncertainties. 

However, as a first step we wish to test whether the method 
works in the distant observer limit with homogeneous sampling and 
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Figure 1. /3 as deduced from mock velocity data (as taken from PSCz high 
resolution data) designed to test various selection effects. The dotted-dashed 
curve shows j3 as calculated from a homogeneously sampled velocity cata- 
logue in the distant-observer-limit. The dashed line shows beta from a ho- 
mogeneously sampled velocity catalogue but with the actual volume cov- 
erage of SEcat, i.e., the distant-observer- limit requirement is relaxed. The 
solid line is /3 as deduced from noise free mock peculiar velocity data with 
the same selection effects as SEcat. The dotted line shows the correct value 
of/3 

noise free data. In this case we have constructed the mock veloc- 
ity data by sampling the PSCz catalogue on every 8-th grid point 
where the velocity is taken to be equal to the z component to mimic 
the distance observer limit. Figure Q shows /3 as a function of the 
smoothing scale (dotted-dashed line) which agrees quite well, es- 
pecially at larger smoothing scales, with the expected value shown 
with the dotted horizontal line. 

Next, the same mock data is used but the peculiar velocity is 
chosen to be the radial velocity, i.e., the distant-observer-limit is re- 
laxed. The points that were chosen for the comparison are restricted 
to the range 30 < r < 60 ft" 1 Mpc from the center of the box (re- 
laxing this restriction alters the results but marginally). The result 
of this test is shown as a dashed line in figureQindicating that the 
recovered (3 is in agreement with its original value. 

The third issue to test is whether the spatial coverage of the 
data set is sufficient, namely, whether the number of SEcat galaxies 
and their actual sky distribution are good enough for a recovery of 
the j3 value? The answer is given by the solid line in figureQclearly 
showing that at small smoothing scales, /3 is underestimated, then 
it increases with the smoothing radius until the correct value is re- 
covered on scales larger than 18 h~ 1 Mpc. Here the positions of the 
mock velocity data are the same as the galaxy positions in the SE- 
cat catalogue but their velocities are taken from the PSCz velocity 
field with no noise addition. 

One might argue that there is no clear convergence of the value 
of p at R s = 24 /i _1 Mpc in fi gureQ therefore, one might need to 
go to larger scales. However, given the size of the current velocity 
data catalouges a larger scale smoothing becomes comparable to 
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Figure 2. as deduced from mock SEcat data where the underlying j3 in 
this catalogues is 0.5 (dotted line). The solid line shows the mean f} as 
obtained from 30 SEcat mock catalogues with the error bars reflecting the 1 - 
a uncertainty around the mean. The dashed line shows the recovered value 
of /3 for the noise free case; clearly, there is an underestimation of /3 from 
the noisy data at small smoothing radii. 



the volume of the data set itself and one has to start to worry about 
sampling issues within the Gaussian kernel itself. As will be shown 
later, the lack of convergence is an issue, however less severe, for 
the real SEcat and PSCz data. 

The final step in our testing is to apply the method to a "full" 
mock SEcat catalogue (with noise and actual sampling). The solid 
line in figure [2] shows the mean value of /3 as a function of the 
smoothing radius as recovered from 30 mock SEcat catalogues with 
the error bars indicating the la scatter about the mean. j3 is clearly 
well reconstructed with large smoothing radii. There is also some 
bias in the mean value of j3 at the smaller smoothing radii with 
respect to the f3 obtained from the noise-free data (dashed line), 
which is probably due to overestimation of the noise variance at 
smaller scales. This is not a big worry as on small scales the PSCz 
catalogue used to produce the velocity data has limited nonlinear 
evolution due to its poor resolution (3.2 /i _1 Mpc) and its velocities 
are purely linear. 

Please note that the error bars shown here are correlated. The 
uncertainty estimates made for this figure, and for the rest of the 
figures in the paper, are based on one of the error bars and not their 
combination. 



3.2 f3 from the real data 

Having tested the method on mock catalogues and demonstrated 
that, on large scales, it gives unbiased results for a SEcat-PSCz 
comparison, we now apply it to the real data. Figure|5|shows the 
measured f3 as a function of the smoothing radius where it has a 
value of 0.6 at smoothing radius of 10 /i _1 Mpc but drops down 
as the smoothing radius increases to ~ 0.45, the curve becomes 
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Figure 3. /3 as deduced from comparison of the real data SEcat data with 
PSCz. The solid line and dashed line reflect results from different assump- 
tions from the density reconstruction. The error bars in reflect the 1-cr un- 
certainty. 

almost flat at R B ^ 18 h~ Mpc. The error bars here are taken from 
the la uncertainties determined from the 30 mock catalogues. 

Branchini et al. (1999) have used two methods to solve for 
the redshift distortion equation (Kaiser 1987) and reconstruct the 
real-space density from the PSCz galaxy redshift distribution, one 
is based on the Yahil et al. (1989) iterative method and the other on 
the Nusser & Davis (1994) spherical harmonic expansion approach. 
In the previous analysis we have used data obtained with the former 
method. However, to examine the robustness of the measured value 
of P we perform the same comparison but with the later method. 
The dashed line in figure [3] shows [3 as a function of smoothing 
radius deduced from the second method which is well within the 
la uncertainty level, albeit being slightly smaller. 

On small smoothing scales the behavior of the curves shown 
in figure |3|is systematically different from those obtained from the 
analysis of the mock catalogues, the former drops with scale while 
the later increases. We attribute this difference to the fact that the 
PSCz catalogue has a limited resolution and its velocity field is 
purely linear. 



4 FUTURE SURVEYS 

4.1 Application to Mock 6dF catalogue 

In the near future, the 6dF Galaxy Survey (Jones et al. 2004) will 
measure the redshifts of around 150000 galaxies, and the peculiar 
velocities of a 15000-member sub-sample, over almost the entire 
southern sky. When complete, it will be the largest redshift survey 
of the nearby universe, reaching out to about z ~ 0.15, and more 
than an order of magnitude larger than any peculiar velocity survey 
to date. Since the two datasets will be obtained from the same sur- 
vey, the galaxy redshift and peculiar velocity catalogues will have 
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the valuable attribute of being subjected to the same selection ef- 
fects. 

Despite the relatively large volume covered by the 6dF galaxy 
peculiar velocity survey the relative nature of the errors in the 
D n — a distance estimation might still diminish the information 
content of the data. To evaluate this effect we apply eq. H3l to mock 
6dF galaxy redshift and peculiar velocity catalogues. In this ex- 
periment the catalogues are constructed from the full nonlinear N- 
body numerical simulation described by Cole et al. (1998), specif- 
ically, the simulation labeled L3S in their paper. The simulation 
assumes a CDM power spectrum of fluctuations with fl m = 0.3, 
A = 0.7, rms fluctuation of the mass contained in spheres of ra- 
dius 8/i _1 Mpc, as = 1.13 and a CDM power spectrum shape 
parameter, F, of 0.25. The simulation box side is 345.6 /i -1 Mpc 
and has 192 s particles. The mock catalogues where produces by 
carving out 6 hemispheres of radius 150/i _1 Mpc of the simula- 
tion box. We obtain the redshift and peculiar velocity catalogues 
with uniform sampling of the galaxies in the simulated hemisphere 
in accordance with the expected sampling of the 6dF survey. The 
real-space distribution is presumed to have negligible errors; but the 
distances in the peculiar velocity catalogues carry errors of 20% of 
their actual values. The input linear bias factor, b is one. 

The star symbols in figure|4|show the average (3 value recov- 
ered from the 12 mock catalogs, as a function of smoothing scale. 
The error bars show the associated variance. The continuous line 
shows the average (3 value recovered from the same 12 mock cat- 
alogues in which no errors have been added to velocities. Clearly, 
the recovered (3 is close to its input value at all smoothing scales. 
If the error level we get is realistic the accuracy with which the 
6dF galaxy survey will recover the (3 parameter (~ 0.05) is indeed 
encouraging. 

The recovery of (3 down to 5 ft _1 Mpc scale is very encourag- 
ing too as it indicates that the Gaussian cell "energy-like" equation 
holds also for the quasilinear regime. Obviously, this point needs to 
be further explored with many simulations and over a wide range 
of point separations. 

4.2 Kinematic and Thermal SZ Clusters 

Inverse Compton scattering of cosmic microwave background 
(CMB) photons off thermal electrons within the hot intra-cluster 
medium of galaxy clusters produce two effects. First, distor- 
tion of the CMB black-body spectrum causing the cluster to ap- 
pear brighter or dimmer at different frequencies and, second, an 
achromatic modification of its surface brightness. These effect 
are known, respectively, as the thermal and kinematic Sunyaev- 
Zeldovich effects (Sunyaev & Zeldovich 1972). The two combined 
with a measure of the cluster temperature give the cluster's radial 
peculiar velocity component to a high degree of accuracy. Cur- 
rent estimates of the measured distance-independent absolute un- 
certainty are as low as 130 kms" 1 (Holder 2004). 

The thermal component of the SZ effect is now routinely mea- 
sured with interferometers and major efforts are underway to sur- 
vey the sky with in the thermal SZ relevant spectral range. Since the 
SZ effect is redshift independent, this kind of survey will provide 
an unbiased catalogue of the massive clusters as far back as their 
formation redshift (see e.g., Carlstrom, Holder & Reese 2002 for a 
review). 

The kinematic SZ is an order of magnitude weaker than 
the thermal component and therefore has been harder to measure 
(Holzapfel et al. 1997; Benson et al. 2003). In the future however, 
the kinematic SZ effect will be measurable and together with the 
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Figure 4. /3 as deduced from mock 6df catalogue as a function of smooth- 
ing radius. The stars show the mean results from 12 independent data sets 
and the error bars show the scatter about it. The solid line shows the result 
obtained from the same analysis performed using a set of noise free-mock 
data. The dotted line is the input value of /3. 

thermal SZ component will provide wide angle surveys of galaxy- 
cluster peculiar velocities up to redshift of about 2. Such data will 
probe the evolution of the dark energy and galaxy-cluster bias evo- 
lution and clearly distinguish between various theoretical scenarios 
of cosmological evolution. 

Like the 6dF galaxy survey, the future SZ surveys will probe 
the density and the peculiar velocity of the same region of space 
with the same objects and therefore allow a measurement of fi m 
(through /3). The left hand panel of figure [5] shows the evoution of 
Hf(Q m ) as a function of redshift for different values of Q m for flat 
ACDM universes normalize to the case of Q. m — 0.3 and S1a = 
0.7 case. This figure demonstrates the sensitvity of the cluster SZ 
peculiar motions to the value of f2 m . 

The right hand panel of figure [5] shows the evolution of 
Hf(Q m ) as a function of redshift for various values of the dark 
energy equation of state parameter, w, in a flat universe (Haiman, 
Mohr & Holder 2001). As pointed out by Lahav et al. (1991) the 
evolution of f(fl m ) partially cancels out with the evolution in the 
Hubble parameter. The weak dependence of the evolution on w 
clearly shows that the equation of state is very hard to measure with 
peculiar velocity data. On the other hand however, this insensitiv- 
ity facilitates a very accurate measurement of the clusters biasing 
factor and its evolution as a function of redshift. 

The expected superior quality of the measured peculiar veloc- 
ity of individual clusters is hampered by their sparseness. There- 
fore, it is essential to analyze the data with methods that are stable 
with respect to this feature. The method developed here is a good 
candidate as it is simple, easy to apply, involves no complicated 
inversion schemes and the vast majority of the measured clusters 
will satisfy the distant-observer-limit assumed in the derivation. 
Initial application of the method to realistic mock catalogues shows 
a good success in the recovery of the (3. However, the mock data to 
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Figure 5. The left hand panel shows the evolution of Hf(Q m ) as a function of redshift for flat cosmological ACDM models with f2 m = 0.1, 0.1, . . . , 1, 
relative to the Q m = 0.3 and Q\ = 0.7 case. The right hand panel shows the evolution of the same quantity as a function of redshift for different dark energy 
equations of state assuming a flat ACDM universe with fi m = 0.3 and Ca = 0.7. The lines are for the dark-energy equation of state parameter, w, values of 
{—0.6, —0.7, —0.8, —0.9, —1. — 1.1, —1.2}, with the thick line showing the evolution in the w = —1 case. Although, the magnitude and extremum location 
of the ratio vary with w, the F function itself depends very weekly on the value of w. 



which we applied it were at redshift zero and limited in size. When 
the cluster mass cut-off exceeds 8 x 1O 13 M0 the simulation box is 
left with very small number of clusters. In order to test the applica- 
bility of the method properly one should apply it to very large scale 
simulations that span the redshift range of — 2; a task that will be 
deferred to the future. 



5 DISCUSSION 

In this paper we introduced the Gaussian cell two point "energy- 
like " equation connecting the two point density and peculiar veloc- 
ity correlation functions. The interpretation of this equation is that 
the change in the velocity correlation function is caused by den- 
sity variation coming from scales larger than the scale set by the 
Gaussian smoothing; this analytic cancellation of the small scale 
power is particular to Gaussian kernels. Two practical applications 
of the Gaussian cell two-point energy-like equation have been de- 
veloped here, the first is direct matter power spectrum estimator 
from peculiar velocity data, and the second is /3 measurement from 
comparison of galaxy peculiar velocity and redshift surveys. The 
later application was restricted to the velocity dispersion, i.e., the 
r = 0, case. 

In the r = case, the relation derived here is similar in its 
mathematical form to the Irvine-Layzer cosmic energy equation. 
This is not surprising as each of the two relations reflect some sort 
of energy balance and should, due to dimensionality arguments, be 
homologous. 

Restricting the main formula to the variance case the relation 
could be easily used to estimate the value of (3 from comparison 
between galaxy peculiar velocity and redshift catalogues. In this 
paper we showed that despite their proximity the PSCz galaxy red- 



shift survey and the SEcat galaxy Peculiar velocity data could be 
reliably used to derive the value of /3. The result is consistent with 
that of previous analyses. The variance case has also been shown to 
apply to the 6dF galaxy survey, despite being far from the distant- 
observer-limit. Using mock 6dF catalogues we have demonstrated 
that our method can be successfully used to extract cosmological 
parameters from the real sample. 

In the future, the eminent detectability of the kinematic 
Sunyaev-Zeldovich effect will provide peculiar velocity measure- 
ments for large number of galaxy clusters at redshifts extending 
back to the formation epoch of cluster (z ~ 2). This type of data is 
ideal to explore with the Gaussian cell two-point energy-like equa- 
tion as it satisfy all of the required assumptions and have small 
measurement errors with large spatial coverage. The redshift cov- 
erage of the SZ data will allow an accurate measurement of the evo- 
lution of the clusters biasing factor with redshift. The main hurdle 
these data sets will pose is the limited resolution with which they 
will sample the universe as the comoving rms distance between rich 
galaxy-clusters is of the order of 30 /i _1 Mpc. 

We have also shown that the Gaussian cell two-point energy- 
like equation could be used to estimate the matter power spectrum 
peculiar velocity data in a non-parametric fashion. This is a very 
important application since the current measurements of the mass 
power spectrum from peculiar velocity employs likelihood analy- 
sis with specific models that almost certainly do not properly ac- 
count for the noise contribution (Zaroubi et al. 1997, 2001). A 
non-parametric measurement on the other hand will allow a scale- 
by-scale dissection of the various components contributing to the 
measured power spectrum allowing the isolation of the noise part. 

Finally, given the simplicity of equation Q it is tantalising to 
attempt to extend it to the quasi-linear regime to obtain a nonlinear 
description of the evolution of the two point peculiar velocity corre- 
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lation function similar to the very successful quasi-linear extension 
of its density counterpart (Hamilton et al. 1991). Indeed, figure |4] 
gives an encouraging indication that the equation might hold for the 
quasi-linear regime. This will be further explored in a future work. 
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